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"i Abstract. The stationary accretion shock instability (SASI) plays a central role in 

O modern simulations of the explosion phase of core-collapse supernovae (CCSNe). It 

■4-^ may be key to realizing neutrino powered explosions, and possibly links birth properties 

j^ of pulsars (e.g., kick, spin, and magnetic field) to supernova dynamics. Using high- 
resolution magnetohydrodynamic simulations, we study the development of turbulence, 

^_H and subsequent amplification of magnetic fields in a simplified model of the post-bounce 

^ core-collapse supernova environment. Turbulence develops from secondary instabilities 

22 induced by the SASI. Our simulations suggest that the development of turbulence plays 

l~ an important role for the subsequent evolution of the SASI. The turbulence also acts 

(^ to amplify weak magnetic fields via a small-scale dynamo. 

rn 

o 

^ PACS numbers: 97.60.Bw, 94.05.Lk, 95.30.Qd 
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1. Introduction 

The core-collapse supernova (CCSN) blast is initiated when the central density of a 
collapsing massive star (> SMq) exceeds nuclear density. At this point the repulsive 
short-range nuclear force stiffens the stellar core's equation of state, and the resulting 
core bounce launches a compression wave that turns into a shock when it reaches the 
supersonically collapsing outer core. As the roughly spherical shock propagates radially 
outward through the collapsing layers of the core, it loses energy through neutrino 
emission and dissociation of heavy nuclei, eventually stalling to form an accretion shock 
some 100-200 km from the center of the star. It is expected that the shock is revived 
within a second or so, enabling it to disrupt the star's outer layers and give rise to the 
supernova. 

In the current paradigm, initiated by Bethe & Wilson [3], the shock is revived by 
energy deposition in the gain region below the shock from the intense flux of neutrinos 
streaming out of the forming proto-neutron star (PNS). However, sophisticated 
spherically symmetric supernova models fail to reproduce the explosion jlQl [311 1^ - 
except in the case of the lightest supernova progenitors with 0-Ne-Mg cores [27] . 
Results from multidimensional multiphysics simulations (e.g., [TH [121 1321 HSl [35]) are 
encouraging, but discrepant results from independent research groups, in part due to 
the physics included in the models and the approximations made, remain to be fully 
investigated. It appears clear, however, that a detailed understanding of the explosion 
mechanism will likely involve the confluence of multi-frequency (or multi-frequency 
and multi-angle) neutrino transport, nuclear physics (equation of state and reaction 
kinetics), rotation, and (magneto)hydrodynamic instabilities [M l Wl \ fIU\. 

The stationary accretion shock instability (SASI [1]) is central to modern 
simulations of neutrino-powered supernova explosions [HI [131 [I2l [321 [121 [35]. The 
dominant SASI modes are often characterized in terms of low-order spherical harmonics 
ym ^£ _ ]^^ ^ _ 0, ±1) [5l |7j. In the early phase, the SASI manifests itself 
through global oscillations of the stalled shock, while in the later stages it results in 
post-shock turbulence (driven by SASI-induced downdrafts from the shock) and an 
overall expansion of the shocked cavity. The instability may facilitate the explosion 
by improving the conditions for energy deposition by neutrinos (e.g., [32]). In terms of 
physical fidelity, the most detailed simulations have so far been carried out in two spatial 
dimensions (2D; with axial symmetry imposed), but early simulations in three spatial 
dimensions (3D) seem to confirm the importance of the SASI (e.g., [211 [16]). Aside from 
the explosion mechanism itself, the SASI may also connect supernova dynamics with 
birth properties of neutron stars (e.g., natal kick [U], spin [6], and magnetic field [15]). 

To further explore the development of the SASI in 3D, we present results from high- 
resolution magnetohydrodynamic (MHD) simulations. We adopt the adiabatic model 
used by Blondin et al. [H [6] and study the development and impact of SASI-induced 
turbulence by perturbing the non-rotating, spherically symmetric initial condition. The 
so-called spiral SASI mode {i = l,\m\ = 1) [HI [71 [18] develops favorably in our 



Turbulence and magnetic field amplification from spiral SASI modes 3 

simulations. Vigorous turbulence develops from secondary instabilities in the shear 
layer separating two counterrotating streams induced by the spiral SASI mode. The 
turbulence feeds on the power in the low-order SASI modes, and our simulations 
suggest that the SASI saturates nonlinearly from the development of turbulence (as 
was proposed by Guilet et al. [22] )• The SASI-driven turbulence is essentially non- 
helical (e.g., [HEl]), and results in an efficient small-scale dynamo. The initially weak 
magnetic fields grow exponentially by turbulent stretching until the magnetic energy 
becomes comparable (locally) to the kinetic energy of the turbulent fiows. The potential 
relevance of our simulation results to multiphysics simulations of CCSNe and neutron 
star birth properties is discussed. 

2. Numerical simulations 

To study the development of SASI-driven turbulence and associated magnetic field 
amplification, we solve the ideal MHD equations (e.g., [29]) for mass density, fiuid 
velocity, internal energy density, and magnetic field (p, u, e, and B, respectively). We 
adopt the ideal gas equation of state, which relates the pressure to the internal energy 
density P = (7 — l)e, where 7 is the adiabatic index. The gravitational potential 
around the central compact object (the PNS) is approximated with the point-mass 
formula $ = —GM/r, where G is Newton's constant, M the mass of the central object, 
and r the radial distance from the center. A time-explicit, second-order HLL-type 
finite volume scheme is used for the time-integration of the ideal MHD equations (see 
[ISlIISlin] for details). 

Our numerical simulations are initiated with the adiabatic setup used by Blondin 
et al. and Blondin & Mezzacappa [HE] (see also [IS]), which resembles the early stalled 
shock phase in a core-collapse supernova: a spherical, stationary accretion shock is 
placed at a radius i?sh = 200 km. A supersonic fiow is nearly free-falling towards the 
shock for r > i?sh- Between the shock and the PNS, the fiow settles subsonically in nearly 
hydrostatic equilibrium. Matter is allowed to fiow through an inner boundary placed 
at -RpNs = 40 km. The adiabatic index is set to 7 = 4/3 (similar to the degenerate 
conditions in the collapsing core). The mass of the central object M = 1.2 Mq and 
the accretion rate ahead of the shock M = 0.36 Mq s~^ are held fixed during the 
simulations. A radial (split monopole) magnetic field is superimposed on the initial 
condition: Br = sign(cos'i9) i?o (-Rpns/''^)^, where Bq is the magnetic field strength at 
the surface of the PNS. We initiate the SASI with small (~ 1%) random pressure 
perturbations in the post-shock fiow. 

We have carried out simulations on a three-dimensional Cartesian grid with 
resolution Al/Rp^^s ~ 0.029, and using up to 1280^ cells (A/ is the width of a cell). Here 
we present simulations where we have varied the initial magnetic field: Bq = 1 x 10^° G 
(model BIO or "weak-field model") and Bq = I x 10^^ G (model B13 or "strong-field 
model"). The initial field can be considered weak in both models in the sense that the 
magnetic energy density is small compared to the kinetic and internal energy densities. 
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3. Turbulence from spiral SASI modes 

The spiral SASI mode emerges favorably from our numerical simulations. Figure [T] 
displays snapshots of the fluid entropy and provides a visual impression of the 
development of the spiral SASI mode and turbulence in the strong-fleld model (the 
weak-field model evolves similarly). Early on (t = 600 ms), the SASI manifests itself 
by small-amplitude global oscillations of the still spherical shocked cavity. In terms of 
an expansion in spherical harmonics Y^™, the shock oscillations are often characterized 
as sloshing modes {£ = l,m = 0) or spiral modes {£ = l,m = ±1). However, sloshing 
and spiral modes are related: spiral SASI modes can be viewed as a superposition of 
sloshing modes out of phase, and thus may be a more general outcome of perturbing 
the spherically symmetric initial condition (non-radially) [TJ UHl HH] • The SASI results 
in angular momentum redistribution and counterrotating fiows inside the shock. 

The spiral SASI mode leads to strong velocity shear in the layer separating the 
two counterrotating streams in the post-shock flow. With growing amplitude, the spiral 
SASI mode eventually leads to the formation of a shock triple point (cf. flgure [I]). 
The shock triple point forms when a shock wave from a steepening pressure wave, 
propagating on the inside the accretion shock, connects with the accretion shock itself 
in the vertex separating the counterrotating post-shock flows and the pre-shock accretion 
flow [6]. The shock triple point (in reality a line segment extending across the accretion 
shock surface) has just formed at t = 700 ms and is located in the top portion of the 
middle panel in flgure [TJ It propagates in the counterclockwise direction. Lower-entropy 
material ahead of the triple point (flowing in the clockwise direction) is diverted down 
towards the PNS, while higher-entropy material trails the triple point. 

Vigorous turbulence emerges as a result of the nonlinear evolution of the SASI. 
The shear layer extending from the triple point down towards the PNS is vulnerable to 
the development of secondary fluid instabilities (e.g., the Kelvin- Helmholtz instability) 
and turbulence. In the right panel in flgure [I] {t = 786 ms), the shock triple point has 
completed one and a half revolutions about the PNS, and the flow appears turbulent 
in a large fraction of the shocked volume. (The low-entropy stream towards the PNS is 
supersonic at this time.) The turbulent nature of the post-shock flow can also be seen in 
the vorticity fleld a; = V x it, which exhibits high spatial and temporal intermittency — 
in part due to the formation of vortex tube structures — in the late stages of the evolution 

[13 [ig. 

The kinetic energy below the shock grows rapidly with time during the ramp-up of 
the SASI. In flgure [2| we plot the total kinetic energy (solid line) in Vgh — the volume 
bounded by the surface of the shock and the surface of the PNS. We also plot the kinetic 
energy associated with radial flow eI^^, non-radial flow E^^^, and "turbulent" flow E^^^ 
(dashed, dotted, and dash-dotted lines, respectively; see caption for deflnitions). 

Saturation of the kinetic energy below the shock coincides with the development of a 
signiflcant turbulent component in the post-shock flow. The kinetic energy beneath the 
shock is initially about 2 x 10^^ B, and the flow is essentially radial until t ^ 400 ms. The 
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Triple point 





Figure 1. Select snapshots showing the development of the spiral SASI mode, and 
subsequent development of turbulence below the shock. The polytropic constant 
K = P/p'^ (a proxy for fluid entropy) is displayed for times t = 600 ms, 700 ms, 
and 786 ms (left, middle and right, respectively). Darker regions are associated with 
higher entropy. 



kinetic energy begins to grow during the initial ramp-up phase of the SASI, which starts 
around 400 ms, and is due to the rapidly growing non-radial component. For t = 600 ms, 
ii^iin is still larger than E^^^_^. However, the non-radial kinetic energy exceeds the radial 
kinetic energy when t = 700 ms, and the two components remain similar in magnitude 
for the remainder of the run. The kinetic energy growth slows down considerably 
around t = 786 ms (cf. figure [I]), but it continues to grow throughout the nonlinear 
phase, with variability on a shorter timescale superimposed. The turbulent kinetic 
energy also grows rapidly during the ramp-up of the SASI, and becomes comparable 
(within a factor of two) to the radial kinetic energy. When averaged over the time 
interval extending from 900 ms to 1100 ms, the total kinetic energy below the shock is 
{E^,X.l: = 0-044 B. Similarly, we find {eIXIi = 0.019 B, {E^JUl = 0.025 B, and 
{E^^X.ll = 0.011 B. Thus, about 25% of the post-shock kinetic energy is associated 
with turbulence (this number, of course, depends on the value adopted for k'^ in the 
definition of E^^^, but we consider /c"^ = 0.1 to be somewhat conservative, and the 
turbulent kinetic energy is possibly even larger; cf. figure [s]). Moreover, the turbulent 
rms velocity m^^^^ = ^J2E^~JM^^ is about 5000 km s~^ during the saturated state (Mgh 
is the mass within Vsh)- 

Additional, quantitative insight into the development of turbulence from the SASI 
is gained by inspecting the Fourier spectra of select quantities from the simulations. 
From the Fourier transform of a vector field /, 

J{k) = ^ [ fix) X exp (tk ■ x) dVL, (1) 

where Vl is the volume of the computational box, we obtain the spectral density per 
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Figure 2. Evolution of the kinetic energy below the shock: total kinetic energy 
^kin = I /v P""^ dV, kinetic energy associated with radial flow E'^i^ — ^ Iv P^r ^^ 
and non-radial flow E^^^^ = 5 /y P("^ + '^m) ^V, and kinetic energy due to small-scale 
"turbulent" flows Sj^ = J^J? ekin(fc) dk, where Ckin is the spectral kinetic energy density 
(cf. ^). In computing E"^-^, we have set k'^ = 0.1 km'^ (i.e., A"^ = 27r/fc'^ « 63 km). 
The diamonds on the solid line correspond to times displayed in figure fl] (only the 
earliest three) and figure [s] (1 Bcthc [B] = 10^^ erg.) 
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where the magnitude of the wave vector (wavenumber) is k = \k\, and ^^2^ is a solid 
angle element in Fourier space. When the spectral density is integrated over /c-space, 
the result equals the real space integral of the square of the corresponding real space 
quantity; i.e., J^"!'"' edk = \ fy |/P dV {k^i^ = 2tt/L and fc^^x = 27r/A/ are given by the 
size of the computational box L and the grid spacing A/, respectively). In particular, 
we obtain the spectral kinetic energy density CMn by setting / = y/pu in ([I]). In a 
similar manner, we obtain the spectral specific kinetic energy density 'e^, the spectral 
enstrophy density Co,, and the spectral magnetic energy density e^^^ by setting / in ([I]) 
to u, uj, and B, respectively. (When computing the spectra, we exclude the pre-shock 
flow, which is not of interest here.) 

The Fourier spectra plotted in figure |3] further illustrate the growth of the low-order 
SASI modes, the development of turbulence, and subsequent SASI saturation. The 
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Figure 3. Evolution of Fourier spectra during SASI development in the strong-field 
simulation. Left 'panel: compensated spectral specific kinetic energy density e^ x fc^". 
Right panel: Spectral enstrophy density e^^. In both panels, dashed lines represent 
time averaged spectra (averaged over the time period extending from t = 1000 ms to 
t = 1100 ms). The dotted vertical reference lines denote spatial scales corresponding 
to 300 km, 20 a;, and 10 AL 



compensated spectral specific kinetic energy density (i.e., e„ x fc^/^) is plotted versus 
wavenumber k for select time states in the left panel in figure IS] (Enstrophy spectra e^^ 
for the same time states are plotted in the right panel in figure [s]) The figure shows the 
time evolution of the spectra from the early, linear phase of the SASI {t = 600 ms; cf. 
left panel in figure [I]), the transition to the nonlinear phase, and well into the nonlinear, 
saturated state. 

Our simulations suggest that the SASI saturates in the nonlinear regime due to the 
development of post-shock turbulence. In terms of spherical harmonics, the SASI is often 
characterized by exponentially growing power in low-order modes (e.g., i = l,m = 0) 
[5]. As a result, the shock surface deviates rapidly from its initially spherical shape. The 
obliquity of the shock front relative to the pre-shock accretion fiow causes the non-radial 
post-shock kinetic energy to grow at the expense of thermal energy |1]. This is consistent 
with our simulations. The growth of the low-order SASI mode is clearly seen in the left 
panel in figure p^ (cf. the evolution of the peak around k = 0.03 km~^ from t = 700 ms 
to t = 900 ms). Growth of the spectral specific kinetic energy density for larger k values 
accompanies the growth of the peak, which reaches a maximum around t = 900 ms. 
The peak in the spectral density is significantly reduced by t = 1000 ms, when the 
system has reached an approximate statistically steady state, and the compensated 
spectral density is nearly fiat over a range of scales (i.e., a Kolmogorov-like spectrum 
Cu oc k~^^^ is established for k G [0.03,0.2] km~^). The dashed line represents the time 
averaged spectrum, which shows a hint of the peak around k = 0.03 km~^. In this 
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saturated state, flows associated with low-order SASI-modes are continuously powered 
by pre-shock accretion, but the turbulent kinetic energy cascade via secondary fluid 
instabilities efliciently saps energy from the low-order modes and prevent them from 
growing further. The turbulent energy is eventually converted into heat via viscous 
dissipation on small spatial scales (or into magnetic energy via a turbulent dynamo, 
and then converted into heat via Joule dissipatioEE] see section 111). Dissipation is clearly 
dominant for k > 0.3 km~^, where the spectral slope of e^ begins to steepen. (Also, the 
enstrophy spectrum turns over around k = 0.3 km~^.) 

A complementary perspective on the development of post-shock turbulence is 
gained from the evolution of the enstrophy spectrum in the right panel in flgure [3] 
(Enstrophy is proportional to the square of the vorticity magnitude Q = u'^/2.) Post- 
shock vorticity is generated baroclinically — and further amplifled by vortex stretching — 
in the strong velocity shear layer separating the counterrotating flows induced by 
the spiral SASI mode |T6]. The vorticity spreads quickly throughout the post-shock 
volume, and is characterized by high spatial and temporal intermittency. (Figure p\ 
shows the spatial distribution of |-B|, which is very similar to that of |ci;| [2J.) For 
t = 786 ms, a growing maximum has formed in the enstrophy spectrum, which peaks on 
small spatial scales. The average wavenumber k^ = J^ '"""^ k'e^^ dk/ J^ '"'"' e^^ dk (marked 
with a diamond on the three most evolved time states in the right panel in flgure |3]) 
is about 0.34 km~^ (A(^ = 1t^ jk^^ ^ 18 km). For later times, the maximum in 
the enstrophy spectrum continues to increase until t ~ 1000 ms, while the average 
wavenumber k^^ remains relatively unchanged (it appears to increase slowly with time). 
We suspect that the location of the maximum in the enstrophy spectrum is sensitive to 
numerical resolution and that k^^ will increase with decreasing A/ (as is the case with the 
magnetic energy spectrum [16]) until some for now unspecifled physical (not numerical) 
mechanism sets in and prevents k^^ from increasing further. Nevertheless, in our 3D 
simulations, the spatial scale associated with the low-order SASI modes (~ 200 km) is 
clearly separated from the spatial scale associated with the turbulence (k^ ~ 20 km), 
and turbulence appears to have a singularly destructive effect on the SASI by draining 
energy from its low-order modes and thereby causing it to saturate. 

4. Magnetic field amplification from SASI-driven turbulence 

SASI-driven turbulence results in post-shock magnetic fleld ampliflcation due to an 
efficient small-scale dynamo [ini [IS]- Here we summarize the main characteristics of 
magnetic fleld ampliflcation in our simulations. Astrophysical turbulence dynamos 
are broadly divided into small-scale and large-scale dynamos. Large-scale dynamos 
(e.g., the afi dynamo) produce large scale magnetic flelds (i.e., the magnetic flelds 
exhibit spatial coherence on scales larger than the flows producing them), and require 
a signiflcant amount of net kinetic helicity (it ■ lj; typically present in shear flows) in 
order to operate (e.g., [33l [9]). We flnd that the turbulence induced by the spiral SASI 

X In our simulations, viscous and Joule dissipation are due to the shock-capturing scheme adopted |16j 
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Figure 4. Magnetic field evolution from SASI-driven turbulence. Left panel: total 
magnetic energy E'^ag — y^ Iv ^^ ^^ (black lines) and "turbulent" magnetic energy 



E' 



-h 



^T t^ina} 



(k) dk (fc^ — 0.1 km^^; grey lines) versus time for the weak-field model 
(dashed lines) and the strong-field model (solid lines) . The results from the weak-field 
model have been multiplied by 10^ for easy comparison of the models. Right panel: 
spectral magnetic energy density Cmag at select times during the evolution of the strong- 
field model. The diamonds on the solid black line in the left panel correspond to 
the times (solid lines) displayed in the right panel. The diamond on the three most 
advanced time states in the right panel indicates the mean magnetic wavenumber /Cmag 
(see text for definition). The dashed line in the right panel represents the time averaged 
spectral magnetic energy density (emag(fc))i;os- 



mode is essentially non-helical, and results in a small-scale dynamo (this characteristic 
may, however, change if a rapidly and differentially rotating PNS is included in the 
model). Small-scale dynamos result in magnetic fields with a characteristic spatial scale 
similar to the turbulent flows that produce them. (See for example Brandenburg & 
Subramanian (TUI for a recent review of astrophysical dynamo theory.) 

The time evolution of the total magnetic energy below the shock is shown in the 
left panel in figure |4] (black lines), for the weak-field model (dashed) and the strong- 
field model (solid). Both models exhibit exponential magnetic energy growth, with 
growth time r ~ 60 ms, during the early nonlinear evolution of the SASI (t < 786 ms). 
The rapid magnetic energy growth in the strong-field model effectively stops around 
t = 768 ms, when the magnetic energy density becomes comparable to the kinetic energy 
density {B ~ ^JJI^u) in localized regions below the shock. The magnetic energy in the 
strong-field model is about 5 x 10^"^ B when it saturates. The magnetic energy in the 
weak-field model grows exponentially with a nearly unchanged growth rate throughout 
the run, and increases by almost five orders of magnitude. 

The evolution of the magnetic energy spectrum is similar to that of the enstrophy 
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spectrum (figure pi). In the right panel in figure HI the spectral magnetic energy density 
from the strong-field model is plotted versus wavenumber for select times during the 
simulation (the same times displayed in figure pi) . A maximum in the spectrum develops 
on small spatial scales (around k = 0.2 km~^), which is due to the developing turbulence. 
The average magnetic wavenumber k^^„ = f, """' ke^^„ dkl f, """^ e„,a„ dk (marked with a 
diamond on the three most evolved time states in the right panel in figure 111) is about 
0.32 km~^ (i.e., similar to k^; A^^g = 2Tr/k^^g ~ 20 km) for t = 786 ms, and does not 
change much for later times. The maximum in the magnetic energy spectrum continues 
to increase until t ~ 1000 ms, when the spectrum is very similar to the time averaged 
spectrum (e„ag)Lo;- 

The magnetic fields induced by SASI-driven turbulence are characterized by a 
highly intermittent fiux tube structure, where all the magnetic energy resides on small 
spatial scales {k > 0.1 km~^). From the grey lines in the left panel in figured we see 
that the growth of the magnetic energy below the shock is due to small-scale fields 
(i.e., i^niag ~ -^lag)- The spatial distribution of the magnetic field magnitude, shown 
in figure [5] from a late stage in the weak-field model, confirms the turbulent nature of 
the 5-fields. The net magnetic fiux below the shock is essentially zero. In figure [5} we 
also see that the magnetic fields resides in a volume well below the shock surface. (In 
addition to exhibiting similar spectral evolution, the spatial distribution of the magnetic 
field and the vorticity field are very similar in our simulations [15].) Moreover, the high 
intermittency (quantified by the kurtosis of the probability density function; e.g., |8]) of 
the magnetic and vorticity fields are very similar, and B and u; also tend to be aligned 
or antialigned. The observed similarities between vorticity and magnetic field during 
the fully developed turbulent state are consistent with the predictions by Batchelor |2] . 
Also, the SASI-driven turbulence is similar in many ways to the convectively driven 
(non-helical) MHD turbulence reported by Brandenburg et al. [S] . 

In our simulations, magnetic field amplification is due to stretching of magnetic 
fiux tubes by the turbulent fiows [15]. The magnetic fiux tubes are "frozen" in the fiuid 
(i.e., fiuid elements remain on the same field line), and, since initially adjacent fiuid 
elements separate exponentially in turbulence (e.g., [38]), the magnetic field strength 
grows rapidly in proportion to the stretching (i.e., the increased relative separation of 
fiuid elements on the field line). At the same time, the fiux tubes undergo a decrease in 
the scale perpendicular to the stretching. In general, the field amplification continues 
until (i) the field becomes strong enough to resist further stretching by the fiuid, (ii) the 
fiux tube cross section becomes so small that further field amplification is prevented by 
resistive dissipation, or, as is the case for our strong-field model, a combination of (i) 
and (ii). 

The turbulent magnetic field amplification is affected by finite numerical resolution. 
(We see from the right panel in figure that the characteristic wavenumber fc^^g lies in 
the diffusive regime of the spectrum.) Both the growth rate and saturation amplitude 
are underestimated by the numerical simulations. We can, however, estimate the 
"true" growth rate indirectly. By adopting a non-ideal electric field with resistivity 
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rj, E = —{u X B) + r/V x B/no, Faraday's law gives an approximate growth rate for 
the magnetic field (e.g., §55 in [29]) 

2" 



1 dB u'^ 



rms 



B dt A^ 



1-R7'^^^ 



■" ' X 



(3) 



where the magnetic Reynolds number i^m = /io u'^^^ A^/r^ has been defined in terms of 
the turbulence scale A^ (equal to A^ or A^^^ ) and the turbulent rms velocity u'f^^. The 
dissipation scale is denoted with A^ (~ A"^ in our simulations). Equation ^ states that 
the magnetic field is amplified if the second term in the square bracket is less than 
unity. Moreover, as long as the magnetic field is weak {u'^^^ is independent of B), the 
amplification is exponential, and the growth rate tends to the inverse turnover time 
for large R^. (The magnetic Reynolds number in the PNS may be as large as 10^^ 
|43j.) With information extracted from the Fourier spectra, we find u'^^^/X'^ ~ 250 s~^, 
which suggests an exponential growth time of a few milliseconds. Since we measure 
an exponential growth time of the magnetic energy of about 60 ms, we conclude that 
the effective magnetic Reynolds number in our simulations is larger than unity (not by 
very much!), but much smaller than physically realistic values. Thus, our simulations are 
underresolved. Simulations with higher spatial resolution can accommodate thinner flux 
tubes, which results in stronger magnetic fields [15]. We also find that the exponential 
growth rate in the simulations increases with increasing grid resolution. 

The presence of amplified magnetic fields does not result in noticeable effects on the 
global shock evolution in our simulations, and the 5-fields are probably unimportant to 
supernova explosion dynamics. This can be understood as a matter of simple energetics. 
Magnetic energy grows at the expense of kinetic energy, and the kinetic energy content 
in the post-shock fiow during vigorous SASI activity (~ 5 x 10^^ B) is not enough for 
magnetic fields to become energetically significant to the explosion (~ 1 B). Moreover, 
the turbulent kinetic energy accessed for magnetic field amplification only amounts to 
a fraction of the total kinetic energy below the shock (figure |2] and figure |4]). These 
observations suggest a rather passive role of the magnetic fields on the overall shock 
dynamics. On the other hand, the sensitivity of magnetic field amplification and 
evolution to numerical resolution prevents us from ruling out completely any effect 
caused by magnetic fields with importance to the explosion dynamics. 

The magnetic fields amplified by SASI-driven turbulence may contribute 
nontrivially to PNS magnetization. By computing the magnetic energy accumulated 
in the volume occupied by the PNS due to the Poynting fiux through the boundary 
surface at r = -RpNs, we estimate that the PNS magnetic field in the strong-field model 
due to this process alone exceeds lO^'' G. The accumulated magnetic energy (> 10^^ erg) 
meets the energy requirements to power the total fiare energy released per SGR and the 
persistent X-ray emission |13] . The weak- field model results in less accumulated energy. 
However, given the issues with finite grid resolution discussed above, it seems plausible 
that small-scale PNS magnetic fields exceeding lO^'* G can form due to the SASI alone, 
independent of the initial magnetic field strength. 
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Figure 5. Spatial distribution of the magnitude of the magnetic field during the late, 
turbulent phase of the weak-field model. Image courtesy of Ross Toedte, ORNL. 



5. Summary and discussion 

Results from three-dimensional MHD simulations of the SASI are presented. The 
simulations are initiated from a configuration that resembles the early stalled shock 
phase in a core-collapse supernova. Although our simulations adopt a simplified physics 
model that excludes neutrino transport, self-gravity, and the PNS itself, the adopted 
spatial resolution is currently inaccessible to state-of-the-art supernova models in three 
spatial dimensions. Our simulations provide valuable insight into hydrodynamic and 
MHD developments, which will likely be relevant to multiphysics simulations of core- 
collapse supernovae. We summarize our simulation results as follows 

1. Vigorous turbulence develops via secondary hydrodynamic instabilities (e.g., the 
Kelvin- Helmholtz instability) in the shear layer separating the counterrotating flows 
induced by the spiral SASI mode, and we find that the specific kinetic energy 
spectrum develops Kolmogorov-like scaling (i.e., e^ oc k~^^^) in the wavenumber 
range separating the post-shock forcing scale (given by the low-order SASI modes) 
and the dissipation scale (set by finite grid resolution in our simulations). 

2. Our simulations suggest that the SASI saturates nonlinearly due to the development 
of post-shock turbulence [22]. In particular, the evolution of the specific kinetic 
energy spectrum (figure Isl) illustrates how the turbulence saps energy from the low- 
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order SASI modes. In the saturated state, a significant fraction of the post-shock 
kinetic energy is due to turbulence (with mJ^^ ~ 5 x 10^ km s~^). 

3. The SASI-driven turbulence results in magnetic field amplification by stretching 
via an efficient small-scale dynamo. The turbulence is essentially non-helical and 
is similar in many ways to convectively driven MHD turbulence (e.g., |8]). The 
magnetic field evolution is sensitive to numerical resolution. In particular, the 
magnetic energy growth rate and saturation amplitude remain uncertain. However, 
estimates using data extracted from our simulations suggest that the magnetic 
energy may grow on a millisecond timescale under more realistic physical conditions 
(i.e., with i?m ^ 1). 

4. The induced magnetic fields have no noticeable effect on the global shock dynamics, 
and this can be understood from simple considerations of the energetics: a 
relatively small fraction of the post-shock kinetic energy is accessed for magnetic 
field amplification. On the other hand, SASI-driven turbulence may contribute 
nontrivially to proto-neutron star magnetization. (Simple estimates suggest PNS 
magnetic fields exceeding lO^'^ G.) Thus, the formation of strongly magnetized 
proto-neutron stars may not necessarily be linked uniquely to magnetorotationally- 
driven supernova explosions (e.g., [30]), but may merely result as a by-product of 
the violent dynamics associated with the explosion of massive stars. 

Given that the SASI may play a central role in facilitating neutrino powered 
explosions, our simulations suggest that turbulence — through its role in the nonlinear 
saturation of the SASI — may play an important role in core-collapse supernova models 
as well. We have intentionally adopted a simplified physics model in order to carry 
out our simulations with high spatial resolution. They therefore represent an initial 
contribution towards a complete understanding of supernova dynamics, and the role of 
turbulence in particular. (Recently, Murphy & Meakin [36j proposed a turbulence model 
for the core-collapse supernova problem.) The role of turbulence in supernova models 
must be further investigated in the context of multiphysics simulations. Fortunately, 
such simulations — in three spatial dimensions — are now becoming available, although 
initially with limited spatial resolution |16]. 

An understanding of the impact of the disparate evolution of turbulent flows in 2D 
and 3D on core-collapse supernova simulations will be important to establish. Current 
state-of-the-art, multiphysics simulations of neutrino-powered explosions pTl [T3l [T2| l32 | 
1121 ES], carried out in two spatial dimensions with axial symmetry imposed, emphasize 
the importance of the SASI in facilitating the explosion. The post-shock flows in these 
simulations are clearly turbulent due to both neutrino-driven convection and the SASI. 
However, it is known that two-dimensional (inviscid and incompressible) turbulence 
supports a so-called inverse cascade, where smaller vortices of equal size coalesce to 
form larger vortices, which leads to increased turbulent energy on spatial scales larger 
than the turbulent driving scale [281 ED]; i-e., opposite to the turbulent cascade in 
3D (cf. figure |3]). Comparisons between 2D and 3D simulations have been carried 
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out in the context of parametrized models |37l [23]. In particular, Hanke et al. [23] 
observed diverging asymptotic behavior with increasing grid resolution in their 2D and 
3D simulations: their 2D models explode more easily with increasing resolution, while 
the opposite is true for their 3D models, and they attribute these differences to the 
different cascading nature of 2D and 3D turbulence. Clearly, a better understanding of 
the role of turbulence in 2D and 3D simulations of the SASI (and how they potentially 
differ) is important to establish. We plan to investigate this issue in detail in a future 
study. 

Neutrino-driven convection develops in the post-shock flow on a timescale that is 
generally shorter than that of the SASI (e.g., [I3j). It will also be important to better 
understand the interplay between the turbulence induced by neutrino-driven convection 
and SASI development — in particular, the influence of pre-existing turbulence on the 
coherence of the low-order SASI modes. 

We cannot completely dismiss turbulence-induced magnetic fields as unimportant 
to the explosion dynamics. With simple considerations of the energetics, we have 
argued above (and elsewhere [151 IB]) that the magnetic fields induced by SASI-driven 
turbulence have little potential to impact the global shock dynamics. The presence of 
strong magnetic fields (model B13 versus model BIO) results in less turbulent kinetic 
energy [TU]. However, given the difficulties associated with finite grid resolution when 
simulating MHD turbulence in large scale simulations such as those presented here, 
one must use caution when extrapolating conclusions based on low Reynolds number 
simulations to systems governed by extremely large Reynolds numbers (e.g., |25j). 
Our simulations merely suggest that MHD turbulence (as opposed to hydrodynamic 
turbulence) may be the more appropriate description of turbulence in core-collapse 
supernovae, and a better understanding of the impact of magnetic fields is clearly 
needed here. Moreover, with respect to magnetic field evolution in CCSNe, we have 
not considered (i) the impact of strong differential rotation near the surface of the PNS 
(which can excite the MRI; e.g., [Il|39]), or (ii) the impact of neutrino-driven convection. 
Studying the confiuence of (i) and (ii) (which alone can excite an afl dynamo [13]) with 
turbulent amplification by the SASI is needed to better understand the role of 5-fields 
in CCSNe. A combination of multiphysics simulations and simulations of simplified 
models (guided by multiphysics simulations) is needed. 
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